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ABSTRACT 

It  is  assumed  that  the  deforrriations  of  a  t  .in  cylindrical 
elastic  shell,  subjected  to  a  uniforra  lateral  pressure,  are  ade- 
quately described  by  the  Donnell  shell  equations.   The  surface 
of  the  shell  is  perforated  by  a  rectangular  cutout.   The  tangen- 
tial raenbrane  dis..laceTiients  and  the  norraal  merabrane  stresses 
vanish  on  the  sim.:ly  supported  edges  of  the  hole  and  the  ends  of 
the  cylinder.   Numerical  solutions  are  obtained  by  replacing  the 
boundary  value  probleri  by  a  finite  difference  approxii.iation  and 
solving  the  resulting  system  of  linear  algebraic  equations  by  an 
iteration  procedure. 

1.  Introduction. 

In  this  paper  vje  study  the  deformations  of  perforated 
cylindrical  elastic  shells  that  are  subjected  to  surface  and 
boundary  forces.   The  presence  of  a  hole  in  the  shell's  vjall 
significantly  alters  the  distribution  and  magnitudes  of  the 
stresses.   A, -proximate  procedures,  which  are  restricted  in  their 
applicability,  have  been  proposed  by  various  authors  [l-6]  to 
determine  these  stresses.   The  paper  by  Savin  [L|.J  also  contains 
a  general  discussion  of  shell  cutout  problems  and  references  to 
previous  work. 

For  a  specific  cutout  problem  we  propose  a  numerical 
method  to  determine  the  stresses  and  displacements  in  the  shell. 
However,  with  suitable  modifications  the  method  can  be  extended 
to  include  a  variety  of  cut-out  shapes  and  other  boundary  and 
loading  conditions. 
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For  the  specific  problem  investigated  the  perforation  is 
assumed  to  be  rectangular,  i.e.  the  ed£,es  of  the  cutout  coincide 
with  generators  and  parallel  circles  of  t^ie  cylinder's  mid- 
surface.   The  shell  is  loaded  by  a  uniform  pressure  applied  nor- 
mal to  the  surface.   On  the  si.iply  supported  boundaries  the  tan- 
gential menorane  displacements  and  the  normal  membrane  stresses 
vanish.   It  is  assuraed  that  the  deformations  are  adequately 
described  by  the  Donnell  thin  shell  theory,  see  e.;^,  [7>5], 

The  method  consists  of  replacing  the  boundary  value 
problem  of  the  Donnell  theory  by  a  finite  difference  approxima- 
tion.  The  resulting  algebraic  equations  are  then  solved  by  a 
special  iteration  i.iethod  involving  tv/o  acceleration  parameters. 
This  procedure  is  described  in  Sec,  3«   The  computational 
methods  required  to  apply  the  method  are  outlined  in  Sec,  I4.. 
In  Sec,  5  the  results  of  the  computation  are  discussed  and  graphs 
of  stresses  and  displacements  are  shown. 

The  ntimerical  results  indicate  that  the  stresses  are  large 
near  the  corner  of  the  rectangular  cutout.   These  high  stresses 
may  have  slowed  the  converti,ence  of  f.e  iterative  procedure.   We 
therefore  conjecture  that  for  problems  v\7ith  smoother  cutouts  the 
method  may  be  more  efficient, 

2.  Formulation.  "  •  •',■•-;    ;   :_••'•-  1 

The  cylindrical  shell  is  of  mean  radius  R,  thickness  t 
and  length  2L.   A  cylindrical  coordinate  system  r,  9,  z  is  em- 
ployed where  r  and  9  are  plane  polar  coordinates  and  z  is  the 
axial  distance  measured  from  one  end  of  the  cylinder. 
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Thus  the  coordinates  are  in  the  ranges  R  -  t/2  J^  r  <  R  +  t/2 , 
-7r<9<7i,  -L<z<L,   The  diraensionless  coordinates  x  and  y 
are  defined  by  the  relations" 

X  5  z/L  ,        |x|  <  1  , 
(2.1) 

y  =  (R/L)9  ,      |y|  <  (tiR/L)  =  j^   , 

The  wall  of  the  shell  is  perforated  by  a  rectangular  hole 
whose  edges  are  on  the  lines  x  =  +  a  and  y  =  +  pj^  where  0  <  a, 
(3  <  1*   The  shell  is  deformed  by  a  iiniforn  pressure  p  applied 
normal  to  the  surface.   The  ends  of  the  shell  and  the  edges  of 
the  hole  are  simply  supported  and  the  normal  membrane  stresses 
and  the  tangential  membrane  displacements  vanish  there.   The 
solution  is  therefore  assumed  to  be  syrrimetric  with  respect  to 
the  planes  x  =  0  and  y  =  0.   Thus  we  consider  the  region  D,  see 
Fig.  1,  v/hich  is  one  quarter  of  the  cylindrical  surface.   Sym- 
metry conditions  must  be  specified  along  the  following  three 
boundary  lines  ofD:  a<x<l,  y  =  0;  0<x<l,  y=/;  x  =  0, 

Pi  <  y   <  i^'  '  • 

The  Donnell  theory  for  cylindrical  shells,  see  e.g.  [7*8], 

can  be  formulated  as  a  boundary  value  problem  for  two  coupled 

fourth  order  partial  differential  equations  in  the  two  dependent 

variables  l.'(z,9)  and  P(z,9),   Here,  \1    is  the  normal  displacement 

of  the  raidsurface  r  =  R  of  the  cylinder  and  EP  is  the  membrane 

"The  coordinates  are  scaled  so  that  the  cylinder  is  of  dimen- 
sionless  length  equal  to  2.   ?Iowever,  there  are  other  scalings, 
say  with  L  in  (2.1)  replaced  by  any  characteristic  dimension  of 
the  shell,  which  yield  formulations  siiailar  to  the  one  given 
below. 
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stress   fiinction  v/here   E   is  Young's  raodulus.      In  dimensionless 
variables   the   differential   equations   of   the  Donnell   theory  are 

(2.2a)  A^w  +  T^§  =1,     /\h  =  T^-f 

where  /\  is  the  biharmonic  operator, 

(2.2b)  £\   -  — r  +  2  x 5  +  — r   . 

On  the  boundary  of  D  we  specify  that 

w  =  — ^  =  f  =  ^C-  =   0    ,      f  or  X  =   1    ,      0   <   y   <  /.  , 

dx'=^  ax^  "    " 

(2.3a)    w  =  ^  =  f  =  ~5  =   0    ,      for  X  =   a    ,      0   <   y   <    p  i', 

ax^  dx^  ■"    ~ 

w  =  ^  =  f  =  —4  =  0    ,      for  0   <  X   <   a    ,      y  =    pi', 
3y2  ay'^  "       " 


chi  _  ^A  _  Jf  _  a^  _  ^ 


(2.3b) 

dw        c\       c^f       d^f 


for  y=0,      a<x_<l,      and 
y=/,      0<x<l, 


-        ^  -         -        3  =  0,        for  x  =  0,      ^k  <  y  <  i. 
3x      dx"^       3x       .-^x"^ 


V 


The  dimensionless  variables  in  (2.2)  and  (2.3)  are 
defined  by 
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w(x,y)  =  aaz,9)  ,   f(x,y)  =    [l2(l-v2)]    "  ^  F{z,Q) 

t 

V  ^3  P   1/2  ,2 

c  =  ■ ^ —  (-h      ,    r=  [i2(i-v2)]         ^    , 

12(l-v^)pL  ^  Rt 

where  v  is  Poisson's  ratio.   The  membrane  stresses  o  ,  o^  and 
c°^  and  the  corresponding  dimensionless  stresses  ^   ,  ^  and  J^^ 

are  given  in  terms  of  the  dimensionless  stress  function 
f(x,y)  by 


I^(x,y)  ^  Ao°(z,9)  =  i^^if-^  ,  q(x,y)  =  A.°(z,9)  =  ^^^  , 


(2.5) 


where 


11°  (^^^,v)  =  Ao°  (z,9)  =    -  ^^^^>y) 
a"-^  =  [I2(l-v^)]    p(L/t)'^  . 


The  quantities  c^,  C/-.  and  c  ,,  ,  which  are  proportional  to  the 

bending  moments,  and  the  dii-iensionless  bending  stresses 

m  ,  m  and  m   are  given  in  terras  of  w(x,y)  by, 
X   y      xy 

m^(x,y)    =  3o!(z,9)    =  slSlliS^  ^    v  '^L^^l      , 

(2.6)  .  : 

m,(x,y)    =  Bo:(z,9)    =  ^!ii%li  .    v  i^-^-ii^      , 

'■» .  — » ... 

...  I 

m^^(x,y)    =    (1-v)3g'    (z,9)    =  £^Jil2LzXL      , 
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where 

t   t  2  • 

B  =  -  ^-{h      . 

12p  ^ 

If  the  displaceraents  of  the  i.iidsurface  of  the  shell  in  the 
z  and  9  directions  are  U(z,9)  and  V(z,©)  respectively,  then  the 
dimensionless  displacements  u(x,y)  and  v(x,y)  are  defined  by 

,    X    ,Rv  U(z.9)      ,    V    /Rv  V(z,9) 
u(x,y)  =    (j-)   ^  '  '    ,      v(x,y)  =  (^)   ^  '  '       • 

C  C 

The  raembrane    stress   displacement   relations   of   the   Donnell   theory 
thus   reduce    to    [7.t6J> 


ro  ^   r  ^^  ^  (dv  ^  ^)-j   v-o  ^  _r  [-41  ^  ^  ^  ^:2hj 
(2.7) 


^^^y    2(i+v)  ^y    ^^ 


The  relations  (2.5) j  (2.6)  and  (2.7)  and  the  conditions 
(2.3)  imply  that  the  required  boundary  and  symmetry  conditions' 
are  satisfied. 

The  complete  forriulation  of  the  boundary  value  problem 
consists  of  the  differential  equations  (2.2)  and  the  boundary 
conditions  (2.3).   The  solutions  depend  only  on  the  parameters 
a,  p,  '[   and  X  •      W'e  obtain  approximate  solutions  for  fixed  a, 
P,  and  ^  and  a  sequence  of  ?' values. 


'The  relations  actually  imply  that  the  membrane  displacements 
tangential  to  each  boundary  are  constants.   However,  these 
constants  are  set  equal  to  zero. 
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3«  Numerical  Methods.. 

To  apply  the  approximate  method  the  region  D  is  covered 
with  a  rectilinear  mesh.   The  mesh  is  of  uniform  spacing  d 
in  the  x  and  y  directions.   The  mesh  is  chosen  so  that  the 
boundary  lines  of  D  coincide  ^^rith  mesh  lines.   There  are  II  +  1 
mesh  lines  in  the  x  direction  and  N  +  1  in  the  y  direction 
so  that, 

(3.1)  6   =  i/n  =  4^1^  . 

The  mesh  points  (x.,y.)  are  defined  by, 

(3.2)  x^  =  i(5  ,  y.  =  j*  ,  i  =  0,1, ...,M;   j  =  0,1, ...W, 

and  the  integers  I  and  J  are  defined  by, 

0 
Xj  =  a  ,  yj  =  p  >:. 

The  mesh  region  Ds,  is  defined  as  the  set  of  points  (x. -y.) 

interior  to  and  on  the  boundary  of  D.   The  mesh  region  d!  is 

o 

obtained  from  Dj(  by  excluding  the  mesh  points:   (x,,,y.), 

j  =  0,...,N;   (^I'yj)*   J  =  0,...,  Ji  (x3L,yj),  i  =  0,...,  I. 

At  each  point  of  D.  the  solutions  w(x.,y.),  f(x,,y.)  of 
(2.2),  (2.3)  are  assumed  to  be  approximated  by  w.  .,  f.  .  which 
satisfy  the  difference  equations, 

(3.3)  A  f •  .  +  ^  L  w.  .  =  0  ,  A  "^  4  +  ^L  f .  .  =  1  ,   (x.,y  .)e  D^ 

Here  the  difference  operator  A,, obtained  by  replacing  deriva- 

o 

tives  in  the  differential  operator  (2.2b)  by  centered  difference 
quotients,  is  given  by 
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L  and  L   are  the  "one-dimensional"  centered  second  difference 
operators  defined,  for  any  mesh  function  Ig.   . i  ,  by 

-2 


Mij  =  (Si.l,j  -  26i,j.  -^  Si.i,j)(6)-    , 


(3.l|b) 


L  g 

y 


-2 


ij  =  (Si,j-^i  -2Si,j  -^  Si,j-in6)-  . 


The  approximating  quantities  jw.  .(  and  If.  .(  must  satisfy  differ- 
ence equivalents  of  the  boundary  conditions  (2.3)  where  all 
derivatives  are  replaced  by  centered  difference  approximations. 
Fictitious  mesh  points  exterior  to  Di    are  introduced  to  apply  the 
difference  equations. 

Thus  a  system  of  2[m(N+1)  -  (I+1)(J+1)J  linear  algebraic 
equations  in  the  2[li(i'I+l)  -  (I+1)(J+1)]  unlcnowns  )W.  .^  and  jf.  •> 

is  obtained,   l.'e  assume  that  as  ^ — >  0  the  solution  of  the 
algebraic  system  converges  to  the  solution  of  {2,2),    (2.3) • 
Hence  for  a  sufficiently  small  S   the  solution  of  the  algebraic 
problem  should  yield  a' close'  approximation  to  the  solution  of 
the  boundary  value  problem. 

For  large  M  and  M,  it  is  impractical  to  solve  the  algebraic 
system  by  a  direct  method,  e.g.  Gaussian  elimination.  Therefore, 
an  iteration  procedure  is  used  to  obtain  approximate  solutions  of 
this  system.   With  w. .   as  an  initial  guess,  for  fixed  values  of 
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a,    p,  /andT,  a  sequence  of  iterates  [w!'^  ,   f."?^]  are  defined 
by  the  recursions 

(3.5) 

Each  iterate  must  also  satisfy  the  difference  equivalents  of 
(2.3).  The  acceleration  parameter  X  should  be  determined  so 
that  the  iterations  {3.S)    converge  as  rapidly  as  possible. 

Therefore,  for  each  n  the  solutions  of  two  algebraic 
systems  of  the  form 

(3.6,     Aj4r  =  *>!;•        K-,'-°; 

are  required  where  the  mesh  function  '.tl^'l    is  either  (zL  \-}^^H 


or 


r-rL^f[^^  +  l]      depending  upon  whether   [s|'^^|  ^^  \^i^-^\ 


°^   Cij   \    '      ^   ■'■^"^  relaxation  procedure,  which  is  a  simpli- 
fication of  the  alternating  direction  method  [9,10]  is  used  to 
obtain  a  crude  approximation  to  the  solution  of  (3.6).   A  single 
horizontal  and  a  single  vertical  sweep  of  the  mesh  is  performed 
instead  of  alternating  back  and  forth  between  horizontal  and  ver- 
tical sweeps,"   Thus  the  a;.proximate  solution  jg.^^"^-'-^<  of  (3.6) 
is  defined  as  tlie  solution  of  the  one  dimensional  difference 
equations  flOl  : 

A  limited  number  of  numerical  experiments  indicate  that  the 
convergence  rate  of  (},S)    essentially  does  not  change  if  several 
sweeps  are  made  through  the  mesh. 
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(3.7)  [(x^,yj)e  d'  , 

Is-^"^    /  and  |g[^"^'''^j  "lust  also  satisfy  the  appropriate 

difference  boundary  conditions.   The  matrices  corresponding  to 
the  one  dimensional  difference  operators  in  (3»7)  are  easily 
inverted  by  factoring.   The  number  p  is  a  second  acceleration 
parameter  and  should  be  determined  so  that  the  iterations  (3.5) 
converge  as  rapidly  as  possible. 


where 


[)■«   Coraputational  Procedures . 

Application  of  the  numerical  procedure  requires  that  the 
mesh  spacing  be  'small'.   An  acceptable  value  of  6    is  determined 
by  comparing  the  solutions  of  a  sequence  of  test  calculations 
with  successively  finer  meshes.   Prom  these  tests  we  conclude 
that  the  6    =  l/2l|.  mesh  yields  sufficiently  accurate  solutions  for 
the  range  of  parameters  considered.   The  numerical  criterion  for 
convergence  is  that  the  iterations  satisfy  the  condition, 

(L..1)    z^^U   max.      ,  ||w(j)-w|-l)Mf(^)-f[-^M|<10-^ 

(x.,yj)  e  D^     "^  "^  '^  ^        ) 

The  value  of  k  that  is  employed  depends  on  the  magnitudes  of 

the  solutions.   For  all  our  calculations  we  used  k  =  8  or  9. 

The  condition,   lim  z    =  0  is  only  necessary  for  convergence. 
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It  is  found,  starting  from  a  given  initial  iterate,  that 
the  number  of  iterations  required  to  satisfy  (J+.l)  depends  on 
the  values  of  n   and  X.   For  fixed  a,  p,  ■t  and  T,  the  optimal 
values  Xrs    and  J^^   are  defined  as  those  values  uhich  i.inimize  t 


ae 


numoer  of  iterations  necessary  to  satisfy  (i^-.l).   To  obtain 
estimates  for  them  a  series  of  test  calculations  are  employed. 
Thus  for  a  fixed  configuration,  i.e.  fixed  a,  p  and^c,    and  usually 
employing  a  coarse  mesh  estimates  of  X   and  p  are  determined  for 
an  increasing  sequence  of  t  values.   These  values  are  used  for 
the  finer  meshes  and  X^  and  p_^  for  other  values  of  2^  are  estima- 
ted by  interpolating  or  extrapolating.   In  Table  I  representative 
values  are  given  of  the  acceleration  parameters  that  v;ere  em- 
ployed. 

TABLE  I 
Acceleration  Parameters  Employed  in  the  Computation 


SHELL   I 

SHELL    II 

r 

X 

P 

X 

f 

31+0 

680 

1020 

1360 

1700 

1.0 
.675 
.30 

.20 
.12 

2[|.6 
9.7 
9.7 
9.7 
9.7 

1.0 
.25 
.10 

10.0 

10.0 

9.6 

9.6 

1 

Numerical  solutions  are  obtained  for  two  configurations, 
which  we  call  Shells  I  and  II,  and  a  sequence  of  values  of 
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~C  =   3I4.O,  6c:c,  102c,  I30O,  17''0.   Shells  I  and  II  are  defined 
respectively  by,  ^  -  ^    *    ^-'  ~  ~22    '    ^~  Ih.     ^^'^  °'  ~     ^^^    ' 
P  =  '/22  ,  y^=   /2I4..   The  initial  iterate  for  a  civen  t'  is 
taken  as  the  'converged'  solution  for  the  nearest  ^T  for  the  sarue 
shell,   i/hen  the  iterations  conver£.;e  the  dimensionless  stresses 
are  coraputed  eraploying  oDvious  finite  difference  equivalents  of 
(2.5)  and  (2.6),   It  is  observed  that  the  number  of  iterations 
necessary  for  convergence  depends  on  the  values  of  the  parameters 
a,  P,  -/  and  ^,   For  example,  with  a  fixed  configuration  the  num- 
ber increases  rapidly  as  6  increases,  and  v;ith  fixed  X^and  ^ 
the  number  increases  as  a  and  p  decrease.   In  all  coses  a  relati- 
vely large  number  of  iterations  is  needed  to  satisfy  (L1..I).   For 
example  using  Shell  I  and  2"=  3i^0,  I4.92  iterations  are  needed,  and 
for  ^  =  1020,  2I4.8I  iterations  are  required.   However,  the  compu- 
ting time  per  iteration  is  relatively  siaall,  e.g.  the  above 
computations  required  a;jproximately  13  and  65  minutes'"  re- 
spectively. 

5»  Presentation  of  Results. 

Representative  graphs  of  the  dimensionless  stresses  and 
normal  displacement  obtained  from  the  numerical  solution  for 
Shell  I  v;ith  '^=   1700  are  given  in  Fig's.  2-8.   In  these  figures 
the  variations  with  the  axial  distance  are  shown  for  a  sequence 
of  y  values.   These  results  indicate  that  near  the  corner  of  the 


'All  computations  employed  v  =  ,3  and  were  performed  on  the 
IBM-7090  computer  at  the  A. E.G.  Computing  and  Applied  Ilathematics 
Center  of  the  Courant  Institute  of  Mathematical  Sciences, 
New  York  University. 
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,,  X  =  a,  y  =  pi;  t 


cutout,  X  =  a,  y  =  p^,  tlie  nvunerical  stresses  are  quite  large. 


Tlie  hole  apparently  afx'ects  the  stress  distribution  everywhere 

p-  0 
The  hoop  membrane  stress  [_ 

.  y 


p-  0 
in  the  shell.   The  hoop  membrane  stress  [_    ^^<i  the  axial  bend- 


ing stress  ^  y.      are  influenced  more  by  the  hole  than  the  other 
stresses.   If  a  =  p  =  0,  tae  stresses  and  displacements  obtained 
from  the  solution  of  tiie  boundary  value  problem  are  independent 
of  y.   These  quantities  are  indicated  with  a  super-bar  and  are 
given  by. 


-•  /  N  _  1    (cos  Px  cor.h  px  +  tan  6  tanh  p  sin  Px  sinhpx ) 

W  ^  X  ^     —    "j     ~  ■ '  '  '  » 

^  T^cos    p  cosh   p(l+tan^p  tanh^p) 


(5.1)  Lx^^)  =Gy(^)-=nxy(^"-)  ^  ^ '  r.°(>^)  =r^(x)  , 


lIm  =^i:Im  -  ^'■"" 


•^  dx 


where 


r2  =  up^  . 

These  results  are  shown  by  the  dotted  curves  in  Pig's.  3,  S,   6 
and  8.   In  Fig.  8,  the  lateral  displacement  w(x)  coincides, 
within  the  resolution  of  the  graph,  v/ith  the  numerical  results 
for  the  line  y  =  0.   These  figures  indicate  that  the  cutout  and 
its  corner  may  increase  the  magnitudes  of  the  stresses  and  dis- 
placements by  large  factors. 
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Although  their  magnitudes  differ,  the  fom  of  the  stresses 
and  displacements  for  other  values  of  X and   for  Shell  II  are 
similar  to  taose  shown  in  Fig's.  2-8,   A  simple  asymptotic  solu- 
tion is  obtained  for  large  2^  by  dividing  the  differential  equa- 
tions by  2^and  letting  T— )  co  while  v;(x,yioo)  and  f  (x,y;oo  ) 
remain  finite.   Integrating  the  resulting  equations  and  applying 
the  boundary  conditions  vie  obtain 

w(x,y}oo  )  =  f  (x,y;oo  )  =  0  . 

To  illustrate  this  result,  graphs  are  shown  in  Pig's.  9-11  of 
typical  stresses  and  displaceiaents  for  Shells  I  and  II  on  tiie 
line  y  =  X  for  an  increasing  sequence  of  T'values.   All  stresses 
and  displacements  seem  to  approach  zero  as  2^ increases.   The 
decrease  in  the  magnitudes  of  the  solution  as  2^ increases  is 
another  possible  cause  for  the  slow  rate  of  convergence  for 
large  2^^, 
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Fig,  1,   The  domain  of  integration,  D,  v;r.ich  is  one  quarter 
of  the  cylindrical  surface. 

Fig.  2,   The  variation  vrith  the  axial  distance  x  of  the  diraen- 
sionlcss  axial  i.iembrane  stress  ^   for  Shell  I  for 
a  sequence  of  y  values.   The  hole  occurs  at  y  =  10/2i\.   . 

Fig.  3»   The  variation  with  x  of  the  dimensionless  hoop  membrane 
stress  T"  for  Shell  I  for  a  sequence  of  y  values.   The 

y 


dotted  curve  is  the  stress  if  no  hole  is  present, 

T-O 
Jii     I'J  X  o  11    y>.     u  J. 

of  y  values; 


Fig,  U-,   The  variation  v;ith  x  of  J|   for  Shell  I  for  a  sequence 


Pig,  5«   The  variation  with  x  of  the  dimensionless  axial  bending 

stress  )   for  Shell  I  for  a  sequence  of  y  values.   For 

^  t 

the  shell  i;ith  no  cutout,  »3  coincides,  within  the 

resolution  of  the  graph,  with  the  x-axis  for  the  range 

0  <  X  <  .833.   For  .833  <  X  <  1  ,   r„  is  represented 

by  the  dotted  curve. 

Fig,  6,  The  variation  with  x  of  ^  for  Shell  I  for  a  sequence 
of  y  values.  For  the  shell  with  no  cut-out,  ^  coin- 
cides,  within  the  resolution  of  the  graph,  with  the  x- 

axis  for  the  range  0  <  x  <  ,875.   For  ,875  <  x  <  1., 

~i  . 

l_^     13  represented  by  the  dotted  curve. 
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Pig.  ?•   The  variation  with  x  of  [^,   for  Shell  I  for  a  sequence 
of  y  values. 

Fig.  8,  The  variation  with  x  of  the  dimensionless  normal  dis- 
placenent  v;  for  Shell  I  for  a  sequence  of  y  values. 

Fig»  9«  The  variation  with  x  of  the  dimensionless  normal  dis- 
placement w  for  Bhells  I  and  II  on  the  line  y  =  ^  for 
a  sequence  of  "z^  values. 


Pig.  10.  The  variation  with  x  of  V°   for  Shells  I  and  II  on  the 
line  y  =  X  for  a  sequence  of  ^  values. 


Pig.  11.  The  variation  with  x  of  V  for  Shells  I  and  II  on  the 
line  y  =  J-   for  a  sequence  of  Rvalues. 
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Attn:     Code  Ii38  (2) 

CO,  Office  of  Hav.  lies, 

Br,  Office,   U9$  Suromer  St, 

Boston  10,  Massac  hisetts  (1) 

CO,  Office  of  Nav,  Res, 

Br.  Office,  J.  Crerar  Library 

86  E,  Randolph  St. 

Chicago  11,  Illinois  (l) 

CO,  Office  of  Nav.  Res. 

Br.  Office,  207  West  2Uth  St. 

New  York  11,  N.Y.  (1) 

CO,  Office  of  Nav.  Res, 

Br.  Office,  1030  E.  Green  St. 

Pasadena,  California         (1) 

CO,  Office  of  Hav.  lies. 

Br,  Office,  Navy  100,  Fleet  P.O. 

New  York,  N.Y.  (25) 

Dir.,  Nav.  Res.  Labs. 

VJashington  25,   D.C. 

Attn:  Tech.  Info.  Officer  (6) 

Code  6200  (1) 

Code  6205  (1) 

Code  6250  (1) 

Code  6260  (1) 

ASTIA,  Arlington  Hall  Station 
Arlington  12,  Virginia        (10) 


Office  of  Techn.  Services 
Dept.  of  Commerce 
Washington  25,  D.C. 


(1) 


Dir.  of  Def.  Res.  and  Engin. 

The  Pentagon,  VJashington  25,  D.C. 

Attn:  Techn.  Library  (1) 

Chief,  Armed  Forces  Special 

V/eapons  Project 

The  Pentagon,  Washington  25,  D.C, 

Attn:  Tecl-in.  Info.  Div.  (2) 

Weapons  bffects  Div.  (1) 

Spec.  Field  Projects  (1) 

Blast  and  Shock  3r.  (1) 

Office  of  Secy,  of  the  Army 

The  Pentagon,   l/ashington  25,  D.C. 

Attn:  Army  Library  (1) 


Chief  of  Staff,  Dept.  of  Army 

\.'ashint^ton  25,  D.C. 

Attn:  Dev.  Br.  (RM)  Div.)  (l) 
Res,  Br.  (RM)  Div)  (1) 
Spec,  lieapons  Br.  (RScD)    (1) 

Office  of  Chief  of  Engineers 
Dept,  of  Army,  V/ashington  25,  D.C. 
Attn:  ENG-HL  Lib,  Br.  Adm, 

Ser.  Div.  (1) 

ENG-VJD  Plan  Div.  Civ.  "..Tcs   (1) 
ENG-EB  Port,  Constr.  3r,, 

Enge  Div,  Mil,  Cons.(l) 
ENG-EA  Struc.  3r.  Eng. 

Div.  iiil.   Constr.        (1) 
ENG-NB  Spec.   Engr.  3r., 

Eng.  RM)  Div.  (1) 


CO,  Engin.  Res.  Dev.  Lab. 
Fort  3elvoir,  Virginia 

Office  of  Chief  of  Ordnance 
Dept  of  Army,  Washington  25,  D.C. 
Attn:  lies,  and  Mato  Jr., 
(Ord.  RM  Div.) 

CO,  Watertovm  Arsenal 
Atertown,  Massachusetts 
Attn:   Lab.  Div. 

CO,  Frankford  Arsenal 
Bridesburg  Station 
Philadelphia  37,  Penna. 
Attn:   Lab.  Div, 

Office  of  Ordnsmce  Research 
2127  Myrtle  Dr.,  Duke  Station 
Durharri,  North  Carolina 
Attn:  Div.  of  Engin,  Sci. 

CO,  Squier  Signal  Lab. 
Fort  Monmouth,  N.J. 
Attn:  Comp.  and  Mat,  3r, 

Chief  of  Naval  Operations 

Dept.  of  Navy,  Washington  25,  D.C. 

Attn:  Op  91 

Op  03EG 

Commsmdant,  Marine  Corps 
Headquarters,  USMC 
Washington  25,  D.C. 

Office  of  Chief  Signal  Officer 
Dept.  of  Army,  V/ashington  25,  D.C, 
Attn:  Engin.  and  Teclin.  Div. 
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Chief,  Buream  of  Ships 

Dept.  of  Navy 

Uashington  2^,   D.C. 

Attn:  Code  106  (1) 

Code  312  (5) 

Code  320  (1) 

Code  370  (1) 

Code  375  (1) 

Code  U20  (1) 

Code  U21  (1) 

Code  U23  (2) 

Code  U25  (1) 

Code  UiiO  (1) 

Code  UU2  (2) 

Code  Iih3  (1) 

Code  525  (1) 

Code  633  (l) 

Cliief,  3uream  of  Aeronautics 
Dept.  of  Navy- 
Washington  25,  D.C, 

Attn:  AE-U  (1) 

AV-3U  (1) 

AD  (1) 

AD-2  (1) 

RS-7  (1) 

RS-8  (1) 

SI  (1) 

AER-126  (1) 

Chief,  .Buream  of  Ordnance 

Dept.  of  "ilavy 

Washington  25,  D.C, 

Attn:  Ad3  (1) 

Re  (1) 

ReS  (1) 

ReU  (1) 

ReS5  (1) 

ReSl  (1) 

Ren  (1) 

Spec.  Proj,  Office,  3ur.  Ord, 
Dept.  of  Navy- 
Washington  25,  D.C. 

Attn:  Missile  8r.  (2) 

Chief,  Bur,  Yards  and  Docks 

Dept,  of  i'Javy 

Washington  25,  D.C. 

Attn:  Code  D-202  (1) 

Code  D-202. 3  (1) 

Code  220  (1) 

Code  D-222  (1) 

Code  D-UIOC  (1) 

Code  D-hhO  (1) 

Code  D-500  (1) 


Commanding  Officer  and  Director 

David  Taylor  Model  3asin 

Viashington  7,   D.C, 

Attn:  Code  lUO  (1) 

Code  600  (1) 

Code  700  (1) 

Code  720  (1) 

Code  725  (1) 

Code  731  (1) 

Code  7U0  (2) 

CO,  UoS.  iJaval  Ordnance  Lab. 
White  Oak,  Maryland 

Attn:  Techn.  Library  (2) 

Techn.  Eval  Dept,         (1) 

Director,  Materials  Lab. 

W.Y.  Naval  Shipyard 

Brooklyn  1,  II.Y.  (1) 

CO,  Portsmouth  Naval  Shipyard 
Portsmouth,  Hew  Hampshire        (2) 

CO,  Mare  Island  Nav,  Shipyard 
Vallejo,  California  (2) 

CO  and  Director 

U.S.  Nav,  Electron,  Lab, 

San  Diego  52,  California         (1) 

Officer-in-Charge 
Nav,  Civ,  Engin,  Res. 
and  ival.  Lab. 
U.S.  Nav.  Constr,  3attal.  Center 
Port  Hueneme,  California         (2) 

Dir.,  Nav,  Air  Experimental  Sta, 
Nav,  Air  Mat,  Center,  Nav,  Base 
Philadelphia  12,  Penna. 
Attn:  Materials  Lab,  (1) 

Structures  Lab.  (1) 

Officer-in-Charge 

Underwater  iixplos.  Res,  Div, 

Norfolk  Haval  Shipyard 

Portsmouth,  Virginia 

Attn:  Dr.  A,H.  Keil  (2) 

CO,  U.S.  Nav,  Proving  Ground 
Dahlgren,  Virginia  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

General  Dynamics  Corp., 

Electr.  Boat  Div. 

Groton,  Connecticut  (1) 
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Supr,  of  Shipbuilding 

USN  and  Nav,  Inspec,  of  Ordnance 

Newport  News  Shipbuilding 

and  Dry  Dock  Co, 
Newport  News,  Virginia         (1) 
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Supr.  of  Shipbuilding 
USrl  and  iJav,  Ins.jec.  of  Ordnance 
Ingalls  Shipbuilding  Corp< 
Pascagoula,  Mississippi 

CO,  U.S.  Nav.  Admin.  Unit 
MIT,  Cambridge  39,  Mass, 

Officer-in-Charge 
Postgrad,  School  for 

n'aval  Officers 
Webb  Inst,  of  Wav.  Arch, 
Crescent  Jeach  Rd, 
Glen  Cove,  L.I.,  W.Y. 

Supt,,  Nav.  Gun  Factory- 
Washington  25,  D,C. 

Comrn,,  Nav,  Ordnance  Test  Sta, 
China  Lake,  California 
Attn:  Physics  Div, 

Mechanics  Div. 

CO,  Nav.  Ordnance  Test  Sta, 
Underwater  Ordnance  Div, 
3202  E.  Foothill  31vd. 
Pasadena  8,  California 
Attn:  Struc,  Div, 

CO  and  Director 

U,S,  Nav.  Engin.  Lxp,   Station 

Annapolis,  Maryland 

Supt,  U.S.  Nav,  Postgrad,  School 
Monterey,  California 

Comm.  Marine  Corps  Schools 

Quantico,  Virginia 

Attn:  Dir,,  MC  Dev,  Center 

Comm  Gen,,  USAF 
Washin^on  25,  D.C, 
Attn:  Res.  and  Dev,  Div, 

CO,   Air  iiaterial  Command 
Wright-Patterson  AF3,  Ohio 
Attn:     MCRLX-B 

Structures  Div, 

CO,  USAF  Inst,  of  Technology 
Wright-Patterson  AFB,  Ohio 
Attnj  Chief,  Appl,  i>iech. 
Group 
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(1) 


(1) 
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(1) 
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(1) 


(1) 
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Director  of  Intelligence 
Headquarters,  USAF  ■ 
VJashington  25,  D.C, 
Attn:   PV  x3r,  (Air  Targ,  Div) 

CO,  A.Fi  Office  3ci,  Research 
Washington  25,  D,Co 
Attn:  Mechanics  Div, 

U,S,  Atomic  Energy  Commission 
l/'ashington  25,  DX, 
Attn:  Dir,  of  Research 

Dir,,  Nat,  Bur,  of  Standards 
Uasliington  25,  D.C, 
Attn:   Div,  of  Mechanics 
Engin,  Mech,  Sect, 
Aircraft  Structures 

Commo,  U.S,  Coast  Guard 

1300  Z,   St,,  iW 

Washington  25,  D.C, 

Attn:   Chief,  Test  and  Dev.  Div, 

U«So  Maritime  Administration 
General  Admin,  Office  31dg, 
i4Ul  G  St.,  iW 
VJashington  25,  D.C. 
Attn:  Chief,  Div,  Prelim, 
Design 

Nat,  Aero,  and  Space  Admin, 
Langley  iiesearch  Jenter 
Langley  Field,  Virginia 
Attn:   Structures  Div, 

Nat,  Aero,  and  Space  Admin, 

1512  H  St.,  NVJ 

VJashington  25,   D.C, 

Attn:   Loads  and  Struc,  Div, 

Director,  Forest  Prod,  Lab, 
Madison,  Wisconsin 

Federal  Aviation  Agency 
Dept,  of  Commerce 
Was'nington  25,  iJ.C, 
Attn:   Chief,  Air  Engin.  Div, 

Chief,  Air  aiid 

Equip,  Div. 

National  Science  Foundation 
1520  He  St.,  NW 
l.ashin^ton,  D.C, 

General  Dynamic  Corp, 
Llectr,  3oat  Div. 
Croton,  Connecticut 


(1) 


(1) 


(2) 


(1) 
(1) 
(1) 


(1) 


(1) 

(2) 

(2) 
(1) 

(1) 
(1) 

(1) 


';j''' 


;'yj..K.i.L!''.  vi''';    'lo  T  .!!:>:':'/! 


\"-U.     »";'3i";  '     "J. I,/    I      ,  'i' 


:f;  ■■.;.■ 


V  loqi-;  •>'•'! 


3i7;r.s.i.-.fii:i!.:    ic  ^"^x^u 


::      ;  '   I  f 


S) 


:;■.;  .VI 


;.   'v 


vo''   ,C,('   (.00 


I  ( *, 


iiO. 


•     »■-    ,< 


.■ui. 


(s) 


.cfe.1  ,fcoi 


(I) 


5f»l.i!30 


i.  ; 


(J.) 


■  vr 


-  u  . 


National  Academy  of  Sciences 

2101  Constitution  Ave,, 

Washington  25,  D.C. 

Attn:  Dir.,  Comm.  on  Ships 

Struc.  Design  (1) 

Exec,  Secy.,  Comm,  on 
Undersea  .;arfare        (1) 

Newport  Mews  Shipbuilding  and 

Dry  Dock  Co, 
Newport  IJews,  Virginia         (1) 

Ingalls  Shipbuilding  Corp, 
Pascagoula,  Mississippi        (1) 

Prof,  Lynn  S,  3eedle 

Fritz  engineering  Lab. 

Lehigh  University 

Bethlehem,  Penna,  (1) 

Prof,  R.L.  Bisplinghoff 
Dept,  of  Aero.  Engineering 
Massachusetts  Inst,  of  Techn, 
Cambridge  39,  iiassachusetts     (1) 

Prof.  H.  1.  Bleich 

Dept,  of  Civ,  Engineering 

Columbia  University 

New  York  27,  Hew  York  (1) 

Prof,  3, A,  3oley 

Dept  of  Civ.  Engineering 

Columbia  University 

New  York  27,  New  York  (1) 

Dr.  John  F,  Brahtz 

Southern  California  Labs 

Stanford  Research  Institute 

820  Mission  St, 

South  Pasadena,  California      (1) 

Dr,  D,0,  Brush 

Struc,  Dept.  53-13 

Lockheed  Aircraft  Corp, 

Missile  Syst.  Div. 

Sunnyvale,  California         (1) 

Prof,  3.  Budiansky. 

Dept,  of  Mech,  Engineering 

School  Appl.  Sciences 

Harvard  University 

Cambridge  38,  Massachusetts     (1) 

Prof.  Herbert  Deresiewicz 

Dept,  of  Civ.  Engineering 

Columbia  University 

632  v.,  125th  St. 

New  York  27,  N,Y,  (1) 


Prof.  D.C,  Drucker,  Chmn, 

Div,  of  Engineering 

Brovm  University 

Providence  12,  Rhode  Island      (1) 

Prof,  John  Duberg 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof,  J.  Ericksen 

Hech.  .engineering  Dept. 

Johns  Hopkins  University 

Baltimore  18,  Maryland  (1) 

Prof,  A.C.  Eringen 

Dept.  of  Aero  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Prof,  \}'^.  Flugge 

Dept,  of  Hech.  Engineering 

Stanford,  California  (1) 

Mr.  M,  Goland,  VP  and  Dir, 

Southv;est  iiesearch  Institute 

8500  Culebra  Rd, 

San  Antonio  6,  Texas  (1) 

Prof,  J,N,  Goodier 

Dept  of  Hech,  Lngineering 

Stanford  University 

Stanford,  California  (1) 

Prof,  L.l;.  Goodman 
Engineering  Experimental  Sta. 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Prof.  M.  Hetenyi 

The  Technical  Institute 

Northwestern  University 

Evanston,  Illinois  (l) 

Prof.  P.G.  Hodge 

Dept.  of  Mechanics 

Illinois  Inst,  of  Technology 

Chicago  16,  Illinois  (l) 

Prof.  N.J.  Hoff,  Head 

Div.  Aeronautical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  VJ.H.  Hoppmann,  II 

Depte  of  Mechanics 

Rensselaer  Polytechnic  Inst. 

Troy,  New  York  (1) 
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Prof.  Bruce  G,  Johnston 

University  Of  liichigan 

Ann  Arbor,  Michigan  (1) 

Prof.  J.  Kempner 
Dept.  of  Aero,  i:ngineering 
and  Appl,  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  Ct. 
Brooklyn  1,  il.Y.  (1) 

Prof.  H.L.  Langhaar 
Dept,  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Prof,  B.J,  Lazan,  Director 
Engineering  ^jqaerinental  5ts. 
University  of  Minnesota 
Minneapolis  lU,  Minnesota      (1) 

Prof.  xL.H.  Lee 

Div.  of  Appl.  Mathematics 

Brown  University 

Providence  12,  Rhode  Island     (1) 

Prof.  George  H,  Lee,  Dir.  of  Res. 

Rensselaer  Polytechnic  Inst. 

Troy,  N.Y.  (1) 

Mr,  S.  Levy 

GE  Electr.  Research  Lab. 

3198  Chestnut  St, 

Philadelphia,  Penna,  (1) 

Prof,  Paul  Lieber 

Geology  Departnent 

University  of  California 

Berkeley  U,  California         (1) 

Prof.  Joseph  I»arin,  Head, 
Dept.  engineering  Mechanics 
College  of  :ngin.  and  Arch, 
Pennsylvania  State  University 
University  Park,  Penna,        (1) 

Prof.  ii.D.  Kindlin 

Dept.  of  Civ,  Engineering 

Columbia  University 

632  v..  125th  3t. 

New  York  27,  New  York 

Prof.  Paul  M.  Naghdi 

Building  T-7 

College  of  Engineering 

University  of  California 

Berkeley  U,  California         (1) 


5  - 
Prof.  VJilliam  A.  Nash 
Dept  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida  (1) 

Prof.  N.M.  Nevmark,  Head 

Dept  of  Civ.  Engineering 

University  of  Illinois 

Urbsina,  Illinois  (1) 

Prof,  E,  Orowan 
Dept.  of  Mech,  Engineering 
Massachusetts  Institute  of  Techn. 
Cambridge  39,  Massachusetts       (1) 

Prof.  Aris  Phillips 

Dept.  of  Civ,  Engineering 

15  Prospect  St, 

Yale  University 

New  Haven,  Connecticut.  (1) 

Prof,  VJ,  Prager,  Chmn, 

Phys,  jci.  Council 

'3ro'.vn  University 

Providence  12,  Rhode  Island       (1) 

Prof.  J.R.M,  Radok 

Dept,  of  Aero  :.ngineering 

and  Apple  Mechanics 
Polytechnic  Inst,  of  3rooklyn 
333  Jay  St. 
Brooklyn  1,  i^I.Y. 

Prof,  E,L,  Reiss 

Inst,  of  Mathematical  Sciences 

New  York  University 

U  Washington  Place 

New  York  3,  i:.Y.  (l) 

Prof,  h.   Reissner 
Dept.  of  Mathematics 
Massachusetts  Inst,  of  Technology 
Cambridge  39,  liassachusetts  (1) 

Prof,  M,A,  5adov:sky 

Dept.  of  Mechainics 

Rennselaer  Polyteclinic  Inst, 

Troy,  New  York  (1) 

Prof,  3,W,  Shaffer 

Deot,  of  Mech,  Engineering 

New  York  University 

University  Heights 

New  York  53,  N.Y.  (l) 

Prof,  J,  Stallmeyer 

Dept,  of  Civ,  i:.ngineering 

University  of  Illinois 

Urban,  Illinois  (1) 
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Prof,  Eli  Sternberg 
Dept.  of  Mechanics 
Brovm  University- 
Providence  12,  Rhode  Island    (1) 

Prof.  T.Y.  Thomas 

Grad,  Inst,  nath  and  iiech, 

Indiana  University 

Bloomington,  Indiana         (1) 

Prof,  S.t^,  Timoshenko 

School  of  -ngineering 

Stanford  University 

Stanford,  California         (1) 

Prof.  A.S,  Velestos 

Dept  of  Civ.  bngineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Dr.  W,  .-enk 

Southwest  iiesearch  Institute 

8$00  Culbera  Rd, 

San  Antonio,  Texas  (1) 

Prof.  Dana  Young 

Yale  University 

New  Haven,  Connecticut        (1) 


Prof.  R.A.  Di  Taranto 
Dept.  of  llech.  ,:;ngineering 
Drexel  Institute 
32nd  and  Chestnut  Streets 
Philadelphia,  Penna. 


(1) 


Mr.  H.K.  Koopman,  Secy, 
Welding  Res.  Council 

Engineering  Foundation 
29  '■'.  39th  St. 
New  York  18,  H.Y.  (2) 

Prof.  ..alter  R.  Daniels 
School  of  Lngin.  and  Arcliit, 
Howard  University- 
Washington  1,  D,C.  (1) 

Comm,,  (Code  753) 

U.S.  Naval  Ordnance  Test  Sta, 

China  Lake  California 

Attn:  Techn,  Library         (1) 

Prof.  J,  ...  Cermak 

Dept.  of  Civ.  Lngineering 

Colorado  State  University 

Fort  Collins,  Colorado        (1) 


6  - 

Prof,  W.J.  Hall 
Dept.  of  Civ.  Lngineering 
University  of  Illinois 
Urban,  Illinois 

Dr.  Hyman  Serb in 
Design  Integration  Dept. 
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